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ABSTRACT 

The phenomenon of thermally sustained pressure oscil- 
lations which frequently occur over liquid helium and 
other liquified gases is modeled and analyzed. The 
model is a distributed system with a linear temperature 
gradient. The interaction between the vapor motion and 
heat transfer is analyzed to derive an expression for the 
time history of the pressure oscillations. 

Cyclic interpretation of the pressure history results 
in a relationship between four parameters which control 
the behavior of the oscillations. The four parameters 
contain information on amplitude of motion, slenderness 
ratio of the tube, characteristic lengths of the tube, 
thermal conductivity, specific heat ratio and viscosity of 
the gas undergoing the oscillations and boundary layer decay. 

The effect of changes within the parameters on the 
theoretical behavior of the oscillating system shows good 
agreement with the behavior of oscillations in previous 
experimental apparatus when sirailiar changes are made. 
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Nomenclature 



The following nomenclature will be used throughout 
the text unless otherwise noted. 



Symbol 

A 

A' 

A w 

*1 

A 2 

a 

C 

C P 

c v 

D 

Dl 

d 

k 

L 

L' 

1 

1 * 



N 1 

N 2 



Definition 

Cross sectional area of tube. 

Surface area for heat transfer. 
Surface area for shear stress. 
Constant coefficient. 

Constant Coefficient. 

Temperature gradient along tube wall. 
Constant . 

Specific heat at constant pressure. 
Specific heat at constant volume. 
Viscous drag coefficient. 

Constant coefficient. 

Constant coefficient. 

Diameter of tube. 

Thermal conductivity. 

Length of tube. 

Reduced length for shear stress. ' 
Length of temperature gradient. 

Length of compressible volume. 

Mass ♦ 

Non-dimensional parameter. 
Non-dimensional parameter. 
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Symbol 

N 3 

n 

P 

P 

*«• 

Q 

R 

T 

Ti 

Th 

Tw 

T 

t 

U 

u 

u o 

V 

v 

X 

x m 

x o 

y 

ym 

r 

s 



Definition 
Non-dimensional parameter. 

Angular frequency. 

Pressure . 

Average pressure. 

Storage vessel pressure. 

Heat transfer rate. 

Gas constant. 

Temperature . 

Cold vapor temperature. 

Warm volume temperature. 

Wall surface temperature. 

Average temperature. 

Time . 

Internal energy. 

Velocity. 

Maximum velocity amplitude. 

Volume . 

Average velocity. 

Direction along tube or flat plate. 

Thullen's maximum position of inertial element. 
Maximum displacement amplitude. 

Height above flat plate. 

Maximum height above flat plate. 

Thermal diffusivity. 

Specific heat ratio. 

Spacial decay rate of momentum boundary layer. 
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Symbol 

s> 

M- 

1 

S 

f 



Definition 

Kinematic viscosity. 

Viscosity. 

Spacial decay rate of thermal boundary layer. 
Density. 

Shear stress on tube wall. 

Non-dimensional parameter. 

Non-dimensional parameter. 



I. INTRODUCTION 



Thermally sustained pressure oscillations are 
often found in the vapor trapped within cryogenic equip- 
ment. Particularly susceptable to this phenomenon are 
the tubes and piping which penetrate storage containers 
of liquified gases. These tubes are characterized by 
being closed at their outside end which is usually at 
room temperature and open at their other end which is in 
communication with the cold vapor of the liquified gas. 
Since it is usually necessary to fill, vent, or measure 
temperature and pressure within storage vessels by means 
of a tube or pipe penetration, construction of cryogenic 
equipment without this configuration is nearly impossible. 

Occasionally the oscillations have been used to good 
advantage for such purposes as vapor mixing or liquid 
level measurement. However, in most ins+ances, the 
presence of the pressure oscillations is undesirable. 

When their pressure amplitudes are sufficiently large, 
the oscillations can make vapor pressure measurements 
erroneous. Since the open end of the tube is in free 
communication with the cold vapor above the liquid, there 
is a heat pumping action to the low temperature region 
which greatly increases the boil-off rate of the liquified 
gas within the container. 

Encounters with the phenomenon in cryogenic apparatus 
is well documented in cryogenic engineering literature. 
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The phenomenon is most readily observable over liquid 
helium. Investigations by Ditmars and Furukawa (3) and 
Thullen (7) have shown that the pressure oscillations 
may occur over liquid hydrogen or liquid nitrogen, although 
their detection is somewhat more difficult than over 
liquid helium. Thullen (7) gives a rather complete survey 
of encounters up to the present time. 

Essentially since the late 1940* s, the literature 
has dealt with the description, detection, and damping of 
the oscillations due to their unwanted nature. Curosiously 
enough, it is Lord Rayleigh (5) that is most often quoted 
in explaining the oscillations. In his work The Theory 
of Sound, Volume II, under the section entitled ’’Maintenance 
by Heat of Aerial Vibrations,” Lord Rayleigh describes 
a form of thermally sustained vibrations encountered by 
glass blowers when blowing a bulb at the end of a long 
narrow tube. Just as found in cryogenic applications, 
the tube was communicating between temperature extremes 
and the tube, itself, certainly must have possessed a 
large temperature gradient along its wall. Qualitatively 
Lord Rayleigh states that, as the gas is moved from the 
cooler regions of the tube to warmer regions, ’’the adjustment 
of temperature takes time (italics), and thus the temp- 
erature of the air deviates from that of the neighboring 
parts of the tube, inclining towards the temperature of 
that part of the tube from (italics) which the air has just 



come. From this it follows that at the phase of greatest 
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condensation heat is received by the air and at the 
phase of greatest rarefaction heat is given up from it, 
and thus there is a tendency to maintain the vibrations.” 
Lord Rayleigh further points out that the bulb and the 
great range of temperature is necessary for the maintenance 
of the vibration. Within cryogenic apparatus construction, 
low conductivity material such as glass, cupro-nickel , 
or stainless steels are used for the tubes. Thus, a 
severe temperature gradient must exist near the cold end 
of the tube; and the warm end of the tube and associated 
valving act as the bulb of Lord Rayleigh's high temperature 
phenomenon. 

Extensive experimentation has been conducted to 
determine the types of materials and geometric configurat- 
ions which will foster and sustain the pressure oscillat- 
ions in a room temperature to cryogenic temperature environ- 
ment. Clement and Gaffney (2), Ditmars and Furukawa (3), 
Bannister (l), and Thullen (7) have in their own ways 
measured frequencies, pressure amplitudes, and devised 
methods of damping the oscillations. From their experiment- 
ions and reports it is evident that the tube diameter, 
length, and magnitude of the longitudinal temperature 
gradient in the tube wall are the most important physical 
characteristics of the apparatus. 

Analytical investigations have yet to lead to a 
satisfactory theoretical explaination of the phenomenon 
associated with cryogenic apparatus. Norton and Mule nhaupt 
(4) have produced a computer program to simulate the 
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pressure oscillations. In their program the continuous 
temperature gradient was approximated by a large number of 
constant temperature regions. The program gave encouraging 
results when compared to experimental measurements of the 
pressure oscillations. Thullen (7) proposed a lumped 
parameter model in which he was able to form two dimension- 
less variables, one containing frequency, friction and 
temperature gradient information, and the other containing 
amplitude and heat transfer information.. From this 
approach the exact behavior of any distributed system 
can not be predicted; but, it does point out the controlling 
variables. Both the approaches made by Norton and Kulen- 
haupt and Thullen have added significantly to the theo- 
retical understanding of thermally sustained pressure 
oscillations in cryogenic equipment. However, due to the 
approaches taken, neither theory was able to capture the 
inherent distributiveness of the system and the continuous 
nature of the temperature gradient along the tube wall. 

It is the intention of this author to define a model 
for the pressure oscillations which encompasses the im- 
portant points of previous experimentation, and analyze 
the model in such a manner as to retain as much distribut- 
iveness of the system and continuity of the temperature 
gradient as possible, while maintaining a manageable 
theoretical discription of the phenomenon. 



II . ANALYSIS 



1 . The Model . 

The model chosen for the theoretical analysis of 
the thermally sustained pressure oscillations closely 
approximates the geometric configuration often described 
in the literature. Common points in these configurations 
are the closed volume at the outside end of the pene- 
trating tube, the longitudinal temperature gradient along 
the tube wall, and the free communication of the cold 
vapor with the open end of the tube. The system for 
analysis is illustrated in Figure 1 . 




Figure 1 



Model For Analysis 
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It is Important that the temperature gradient 
along the tube wall be of sufficient strength to allow 
heating of the vapor upon compression and cooling upon 
expansion. The maximum amplitude of the oscillations 
is limited due to viscous friction. Thus the model of 
real systems has all the elements to foster the growth of, 
and sustain, pressure oscillations. Also, the distributed 
nature of the system and the continuity of the temperature 
gradient have been preserved. 

Within cryogenic apparatus the oscillations are 
spontaneous in origin; and, after a period of time, their 
amplitudes grow to a maximum value. The analysis of this 
model will assume that a sinusoidally varing steady state 
condition has been reached within the system. The re- 
lationships between time, motion, and heat transfer will 
be combined to define the pressure as a function of time, 
physical characteristics of the vapor, and geometric 
configuration. 

2. Derivation of the Time and Pressure Relationship. 

When the model is considered as an open system, 
the oscillations of the vapor within the tube cause 
energy to be added or lost to the system by virtue of 
mass entering or exiting the open end of the tube, and 
by heat transfered to the vapor because of the temperature 
gradient along the tube wall. The first law of thermo- 
dynamics for an open system applies to the control volume 
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within the closed tube: 

dU _ . * dm , * , . s 

dt ~ C P T ldt + ^ ' 1 ' 

The time rate of change of mass entering the system 
across the open end of the tube is given by: 

ff =/Av(t) (2) 

where v(t) is the average velocity of the fluid particles 
at the end of the tube at a given time t. 

The perfect gas relations and the definition of 
Internal energy give the rate of change of pressure in 
terms of the rate of change of internal energy: 



p mRT 


(3) 


V 


U = mc v T 


(4) 


p = JL u 


(5) 


C v V 


O 

II 

P* 


(6) 



where C = —— 
c v V 

Since both pressure and internal energy are functions 
of time : 

dP _ r dU 
dt “ dt 



( 7 ) 
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When equation (7) is substituted into equation (l), 
the rate of change of pressure is given as a function of 
time : 



Density is not a variable of primary interest; and, 
over a cycle of the oscillations, its variation is small. 
Therefore, the following simplifing assumption and sub- 
stitution is made in equation (8): 




where P is the average pressure at an average temperature, 
T, during a cycle; 



When v(t) and Q, ( t ) are known, equation (10) can be 
integrated to give the pressure as a function of time. 

Since the usual tube diameters are large in comparison 
to the thermal and momentum boundary layer thicknesses, 
the derivation of the velocity term, v(t), and the heat 
transfer term, Q,(t), in equation (10) will be carried out 
in a rectangular co-ordinate system. Appendix A shows the 
range of tube diameters for which this approach will be 



— = C^c p A^T 1 v(t) + Q(t)] 



( 8 ) 



dP = c 
dt 




( 10 ) 



valid 
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3 . Solution of the Velocity Term , v ( t ) . 

Schlichting (6) presents the boundary layer solution 
for the flow of a fluid near an oscillating flat plate. 

For the purpose of this analysis, Schlichting' s solution 
will be modified to fit the boundary condition that the 
fluid is sinusoidally oscillating over a stationary flat 
plate. Accordingly, the velocity distribution in the 
fluid above the flat plate as a function of height above 
the plate, y, and time, t, is given by: 

u(y,t) = u 0 e"^ycos (nt - <fy) - u 0 cos nt ( 1 1 ) 

where £ =~\f is the spacial decay rate of the momentum 
boundary layer. 

The average velocity of the fluid entering the tube 

is obtained by assuming that the average velocity of the 

fluid over the flat plate to a given distance is analogous 

to that entering the tube: 

r ym 

vftl = Jou(y,t)dy (12) 

y m 



where y m corresponds to the tube radius, d/2, and is 
sufficiently large so that in equation (12) e”^ m <§d1. 

When the average velocity is determined in this manner, 
the final solution of the velocity term is: 



v(t) = u 0 



— .1 _ - l)cos nt + 1 sin nt 

2<*y m 



( 13 ) 
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4. Solution of the Heat Transfer Term , Q. 

The model assumes that the only heat transfer to and 
from the vapor Is in the area of the tube where the 
temperature gradient exists. The rest of the volume of 
the tube is assumed to be adiabatic with respect to the 
vapor. As with the velocity term derivation, the flat 
plate approach is assumed to be applicable. 

To begin, an energy balance is made on a control 
volume of unit depth positioned above the flat plate. 
There is a vertical heat flux through the control volume 
and a horizontal enthalphy flux. Figure 2 illustrates 
the energy balance. 




77777777777777777777777 
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Figure 2. Energy Balance For Heat Transfer Derivation 



The energy balance on the control volume yields the 



1 1 



following differential equation: 




(14) 



Equation (14) may be rearranged into a diffusion equation 
with a velocity type forcing function: 



When equation (15) has been solved for the temperature 
distribution in the fluid above the flat plate, the 
temperature gradient at the wall can be found by differ- 
entiating the solution with respect to y and setting y 
equal to zero. This will enable the heat transfer term 
to be realized as a function of time. 

The following assumptions are made in reguard to the 
solution of equation ( 1 6 ) : 

a. A linear temperature gradient exists in the wall 



b. The same linear temperature gradient as in the 
wall exists within the gas so that the form 
of the temperature distribution in the gas above 
a point on the plate is given by: 



1 dT _ ^ tl. - u(y,t)lT 

2^ c>t 3 y 2 & x 



(15) 



where d* = -J- and cf - 



k 




so that the wall surface temperature is given by 



T w = ax + E 



( 16 ) 



T(y, t) = f (y , t) + ax 



(17) 
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c. There is neglible heat conduction in the gas 
in the x-direction. 

d. The velocity distribution is given by equation 
(ll) and the average velocity is given by 
equation (13). 

Equation ( 1 7 ) must also fit the boundary condition 
that at the wall the wall surface temperature and the 
gas temperature must be equal. 

When equation (17) is substituted into equation (15), 
the differential equation which determines the temperature 
distribution above the plate is: 



1 *1 
7* at 




au(y,t) 



( 18 ) 



where the function f(y,t) is subject to the boundary 
condition that at y = 0, f(0,t) = E. 

The solution equation ( 1 S ) is of the form: 

f(y,t) = A^sin (nt - <^y) + A 2 sin (nt - iy) 

+ B.sin nt + B^cos nt 

1 2 

+ e"^ D^sin(nt - ty) + D^cos^t - f y) 

+ constant (19) 



Appendix B shows the determination of constants associated 
with equation (19). 

The complete solution to equation (17) giving the 
temperature distribution in the gas above the flat plate is: 
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T(y,t) = E + ax + u 0 a e - ^ sin(nt - <Sy) 

2 c<<$ i - n 

7 * 

+ u 0 a3^sln nt 
n 

- 2u Q ac<(£ Z ?/ / n e"^ysin(nt - fj) (20) 

2^cT 2 - n 

where f -yf ~ n , spacial decay rate of thermal boundary layer 

2* I 4 . 

S =~\j n , spacial decay rate of momentum boundary 
v 2 *v layer. 

The heat transfered from the plate to the oscillating 
vapor is given by: 

k = -k^T| A* ' (21) 

c>yjy =o 

When equation (20) is differentiated with respect to 
y and evaluated at y = 0 and the result substituted into 
equation (21 ), the heat transfer term, Q,, is given by: 

Q, = - kaA 1 Upcf \.2<*S fd 1 - 1 
2 C {<$*- n [ n 



5. Fina l Solution of the Pressure-Time Relationship . 

With the results of sections 3. and 4. combined 
with the differential equation of section 2., the final 
derivation of the gas pressure as a function of time can 
be determined. 

Equation (10) of section 2. is: 



f 



in nt + cos nt 



( 22 ) 
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dP _ c 


CpPA 


3.1 


dt 


1 R 


T 



v(t) + Q,(t) 



Equation ( 1 3 ) of section 3. is: 



v ( t ) = u 0 



2<fy 



- 1 



m 



cos nt + _i_ sin nt 



Equation (22) of section 4. is: 



1" _ I 


sin nt + cos nt 


L n 


Jl J 



When equations ( 1 3 ) and (22) are substituted into (10) 
and equation (10) is integrated with respect to time, the 
gas pressure as a function of time is: 



P(t) = C, 



(s#| 


[a) 




1 R / 


w 


n JL' 


f kaA 1 & 




Hcj 


1 2«'<f a - 


n 


»1 



1 _ ij sin nt - 1 cos nt 

|2<fy ra I 2<5y ffi 



n 



+ P e 



-cos nt + sin nt 



(23) 



} 



where P«> is the constant of integration or pressure at 
time t = 0”. 

With the following definitions the compressible 
volume pressure given by equation (24) can be made non- 
dimensional and more manageable for practical applications: 
Let 



N 0 = —L 

° 2 a 



Ym 

N 1 = 2 JSfz 4 



n 



N 2 = 2 

n 
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T = h 

T 

H = 

R 

B = kaA'J 
n 

Upon the substitution of the above definitions, 
equation (23) becomes: 



Pit) 

X BCu 0 /ri 




HT* 

B 



N 



o 




N4 - 1 
2 " 




sin nt 



cos nt 



In a likewise manner the average velocity, v(t), 
be made non-dimensional: 



v(t) 

u 0 



(N 0 - i) cos nt + N 0 sin nt 



( 24 ) 

can 



( 25 ) 



III. DISCUSSION 



With the pressure history of the warm volume and 
the motion of the mass taking part in the oscillations 
now known, a dynamical approach to the system will be 
used to derive the relationship between the parameters 
which control the oscillations. The model of the system 
for this Investigation will be similar to that used by 
Thullen (7) and is shown in Figure 3. 



— c*| d |«s — 





Shear Stress 
on Tube Wall 



Figure 3* Model For Cyclic Interpretation 
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Newton's law of motion applies to the inertial 
element of the model: 

mx + Dx - A ( t ) - P e ( t )j = 0 



(26) 



The transient solution to equation (?6) is not of 
interest and can be avoided by the cyclic integration 
of equation (26) with respect to displacement. The 
net work done in accelerating the inertial element through 
a cycle is zero: 

^ mxdx =0 (27) 

The two remaining terms of equation (26) are now 
equated and point out that the cyclic work done on the 
volume to the left of the inertial mass is equal to the 
cyclic work done against friction: 



^A^P(t) - P e (t)j dx =^Dxdx 



( 28 ) 



The identity dx = dt = v(t)dt is substituted 

dt 

into equation (28) so that the cyclic integration can 
be performed with respect to time instead of displacement. 

A ^P ( t ) - P e (t)jv(t)dt =^Dv 2 (t)dt (29) 



/ 



If the exit pressure is constant, P e (t) = Pea, the 
following eneauality must be satisfied for the oscillations 
to be self sustaining: 



^A^P(t) - Pcoj v(t)dt = ^>Dv 2 (t)dt 



(30) 
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In substituting equations (24) and (25) into (30) 
and performing the cyclic integration, the following 
relationship is established: 



ABC 

nD 



N i - 1 

N 2 - 1 



£(K 0 - I ) 2 + NS 



(31) 



Thullen's (7) non-dimensional parameter containing 
frequency, friction, and temperature gradient information 
is as follows : 



N t = Zcfi 2 . b _i_ (32) 

s D coTT 

Upon substitution of the following definitions 
of the variables in equatiion (32), the Thullen parameter 
becomes : 



P 0 = /RTx 



A 



= TTd2 



S = TTd^l * 

T" 



b 




co = n 




( 33 ) 
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Upon sutsitituion of the actual variables for B and 

C In the non-dimensional grouping ABC in equation (31 ) and 

nD 

making the following definitions, this non-dimensional 
group is seen to contain Thullen's parameter plus additional 
information : 



N* = ABC = AkaA'J V R 
nD n^D c v V 

A = TTd 2 

T~ 

A' = TTdl 

V =7fd 2 1’ 

T* 

c f =’\f'rT = ~\fnf = j>-\ T n 

2* 2yU. 2/t/> 



4,-rrk 

c v r.^ dyb/t/^ 

N* = N T F 4 a rk 1 

[5^la yWj 



N* = 






_i 

Dn 



(34) 



(35) 

(36) 



In order to carryout the cyclic integration of 
equation (2?) with a variable exit pressure, Pe(t), 
further definition .of the exit pressure term and the 
damping coefficient term, D, are required. 

First, the exit pressure is examined: 



Z 



pit) 






( t ) 



p» + v(t ) 

Outflow: P e (t) = P«o 
Inflow: P G (t) = Pco - i/v 2 ( t) 
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Q> AP e (t)v(t)dt = /AP e (t)v(t)dt 

J A - 



(37) 



nflow g- cycle 

+ f AP e (t)v(t)dt 
-'out flow -§■ cycle 

The outflov/ half-cycle contributes nothing to the 
cyclic integration. In the inflow half-cycle the sign 
of the velocity is negative. Therefore, the cyclic 
integration of equation (29) reduces to: 

<£ A Fp ( t ) - P ftJ lv(t)dt - f-|Aj>v^(t)dt = <^Dv2(t)dt (38) 

^ -'inflow f cycle 

Second, the damping coefficient term can be determined 
through a pipe flow analogy involving shear stress: 



E 



3 



F = (P 1 - P 2 ) A =TwA w 



7T W - u sLU 
^ c>y 



y = o 



The rate of change of velocity with respect to y 
is found from the velocity equation ( 1 1 ) 



<Du 

dy 



Q = -(/uc^sin nt +cos ntj 



(39) 



The damping coefficient is now defined as: 

D = -F = A w ^(u 0 [sin nt + cos nt] 
v(t) ' “(t ) 



(40) 



The right half of equation (38) now requires the 
cyclic integration of the following term: 



X + cos nt ] v 2 ( t ) at 



(41) 
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With the above substitutions for the exit pressure 
and damping coefficient terms and subsequent cyclic 
integration of equation (38), the relationship between 
the parameters which contribute to the sustaining of the 
pressure oscillations can be established: 



ABCugTT 

n2 



N i - V 
N 2 - 1, 



“ [< N o “ !) 2 + *o] 



3/2 



1/2 



irA w _k_u 0 [(k 0 - l) 2 + N 2 ] (42) 



Upon substitution of the following definitions 
and previous values of B and C, equation (42) takes 
the form: 




fe _ z _ L | 

1*2 - 1 / 



= _hL_ 
l oty 



-i) 




f = (*c - 








u 0 = nx 0 

A = trdL* , where L' is the reduced length (L - ■§■1*) 
of the tube since the shear stress is 
not uniform over the entire length of 
the tube of Figure 3. 




( 43 ) 
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Figure 4. shows the non-dimensional displacement 

amplitude, 1 4V— olotted against the driving force 

3 U 2^ 0 

oarameter, — l) with lines of constant Uf . 

l'L>n 2 / 

The graphical representation of equation (43) gives 

indication of a given system's performance and how that 

performance will change if the various parameters are 

altered. The manner of change suggested follows closely 

that which has so often been described in literature on 

thermally sustained pressure oscillations. Due to the 

lower limit on tube diameter imposed by the theoretical 

approach taken, the lines of constant JJ in Figure 4. 

are in aggrement with the theory for Iff = 1 . Values of 

9 

\jf = 1 are shown in dashed lines for completeness. 

Of particular interest is the linear relationship 
between the displacement amplitude and the slenderness 
ratio, L'/d. Bannister (l) noted this in his experimentation 
with oscillations over liquid helium. He concluded that 
the pressure amplitudes are directly proportional to the 
slenderness ratio of the gas column. The Bannister data 
was used to calculate values of \jf and the driving force 
parameter. The calculated operating points are plotted 
on Figure 4. Appendix C tabulates the constants and data 
used. Figure 5. shows the calculated amplitude parameter 
plotted against the measured maximum pressure amplitude. 

The actual slenderness ratio of Bannister's tubes is also 
plotted against the maximum pressure amplitude to show 
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Figure 4. Study of Oscillation Amplitude 
Parameter verus Driving Force 
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the linearity comparison. 
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Figure 5. Slenderness Ratio Comparison 

Thullen (7) also found a non-dimensional displacement 
amplitude parameter in his analysis: x m /21. The char- 
acteristic length, 1, was the distance between two ideal 
heat exchangers in his model; and, it was viewed as the 
distance traveled by the gas after a flow reversal before 
significant heat transfer begins. However, the length 
was not well defined for a distributed system. The 
results of the present analysis indicate that Thullen' s 
"1” might have the following value for a distributed system: 

i = 5 j'VET 

d n 
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IV. CONCLUSIONS 



A greater understanding of thermally sustained 
pressure oscillations in liquified gas apparatus is 
possible through a study of vapor motion and heat transfer 
interaction. From such a study a time dependent expression 
for the pressure history of a vapor in a tube can be 
derived. With cyclic interpretation the relationship 
among the control ing parameters of the system can be 
established . 

When changes are made in the parameters, the behavior 
of the theoretical relationships appear to follow closely 
that which occurs when slmiliar changes are made in 
experimental apparatus. 

The rectangular co-ordinate system from which the 
motion and heat transfer interaction was derived places 
a lower bound on the diameter of the tube for which the 
theory applies. Therefore, an analysis with the same 
general approach to the motion and heat transfer, but in 
cylinderical co-ordinates, would be of value. It is not 
anticipated that this would alter the basic relationship 
between parameters controling the oscillations; however, 
the cylinderical system might provide more exact numerical 
constants. Appendix D suggests an approach to the analysis 
of the phenomenon of thermally sustained pressure oscillations 
in a cylinderical co-ordinate system. 



APPENDIX 
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APPENDIX A 



In order to show that the flat plate analogy is 
valid for the results desired, it is sufficient to show 
that the spacial decay rate is small in comparison to 
the tube radius. 

r = 

o 2 S> 

Let <5 y = y/x 

n 

5 x 10" 4 ft ? /sec 
50 rad/sec 

then A = 0.0535 in 




Hence, tubes with radii>0.125 in. will be in the 
range for which the theoretical approach is valid. 
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APPENDIX B 



The determination of the constants associated with 
equation (19) which is the solution to the differential 
equation ( 1 8 ) is as follows: 



Equation ( 1 8 ) : » ^f 

T <*t 




au(y,t) 



where u(y,t) 



u o e 



~ J y C os (nt - c5y ) 



- u 0 cos nt 



and f(y,t) must satisfy the "boundary condition: f(0,t) = E. 
Equation ( 1 9) : 

f(y,t) = e'^^A^ sin ( nt - cfy) + A 2 cos(nt - cfy)] 

+ B^In nt + BgCos nt 

+ e"^y^D^ sin(nt - ? y) + DgCos(nt - f y)j 

+ constant. 

Thus the form of the function f is: f(y,t) = fj + fg + Constant 
Seek the solution of equation (18) by f ^ : 

= e ^|Ajncos(nt - <$ y) - A 2 nsin(nt. - <fy)j 
t Buncos nt + Bgnsin nt. 

= 2 6 Z e~^ ^-A 2 sin(nt - Sy) + AjC03(nt - £y)J 

Upon substitution of these partial derivatives into 
equation ( 1 8 ) and with like coefficients of the sine and 
cosine terms equated: 



A < -■ una 

1 (2of cT 2 - n) 
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A 2 (2c«*cf 2 - -) = 0, A 2 = 0 
_ 

Bi = uoa-rf 
B2 = 0 

f i (y,t ) = Ho5®_f^sin(nt - <T y ) + uoaloos nt 

2*S a - R n 

3" 

However, fj, alone, will not fit the boundary condition. 
Hence, an additional function must be added to f^ and 
that function must be a solution to the homogeneous part 
of equation ( 1 8 ) : 1 ^f _ -t & f 

7t at - d y 2 

f 2 (y» t ) = e ~^ y [ D isin(nt -fy) + DgCos^t -fy)] 

+ constant 

= ne"^y J-D 2 sin(nt - fy) + Djcosfnt - fy)] 

5 ^- 2 = 2 f%~ y[-D 2 sin(nt - fy) + D 1 cos(nt - fy)] 

Upon the substitution of the partial derivatives 
into the homogenous differential equation and with like 
coefficients equated, f 2 becomes: 

D 2 (2*f 2 -^) = 0 Dj(2 = 0 

If Di / 0 and D P ^ 0, then (2<^f 2 - £) = 0 

c 7r 

Therefore, , 

f = V — 

( V 2°t& 

V/ith the boundary condition now applied to f(y,t), 
the cceffients D^, Dg and the ’’constant" of f 2 are determined: 
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constant 

Therefore , 

constant - E 



2u o&w n 



2 **<r- A 



2^ J 



+ DA sin nt + D 2 cos nt = 



d 2 = 0 



'1 ~- 



_ 2u a»cj ft 



2*cf*- B, 



f 2(y*t) = e-ry 



-2u a 






2< 3 f6 s - £ / 

a" 



sin nt + E 



APPENDIX C 



The following constants for helium gas (-200°F) 
were used in the Bannister data fit of equation (43) : 

of = 5.7 x 10 “* ft 2 /sec 
yU. = 8.3 X 10- 6 lbm/ft-sec 
= 4.0 x 10"^ ft 2 /sec 
k = 2.2 x 10-5 bTU/ sec-f t-°F 
= 1.67 



Data from Bannister’s exoerimentation : 





n (cps) 


d (cm ) 


L (cm) 


P (atm x 


10 3 ) AT (°F) 


X 


33.0 


.287 


149 


56 


480 


o 


33.4 


.482 


149 


28 


480 


□ 


33.4 


.605 


149 


18 


480 


A 


22.7 


1 .240 


213 


12 


480 



Estimated from apparatus geometric configuration: 





L' (cm) 


1' 


X 


105 


89 


o 


105 


89 


a 


105 


89 


A 


105 


89 
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APPENDIX D 



The cylinderical co-ordinate analysis will proceed 
with analogous steps to those performed in rectangular 
co-ordinates. Figure 5. illustrates the system for analysis. 
Schlichting* s solution for the velocity distribution of 
an oscillating flow through a pipe is: 



Navier-Stokes Equation: 



du _ 1 SP , .v/6^U , 1 chu 

it - t »x 2 + ?s?> 

Assume that the pressure gradient is caused by a harmonically 
moving piston: 



1 ^P _ v „4- 

f> bx " KcoS nt 

Use complex notation and attribute physical significance 
only to the real part: 

- j|| = Ke int 

Assume the velocity function has the form: 



u(r,t) = f-(r )e 



int 




Figure 6. Cylinderical Co-ordinate Model 



32 



33 



The solution to the Navier-Stokes equation which 
gives the velocity distribution within the pipe is: 



u(r,t) = - i^e iri ' t 

where J 

-y) 



<To(r <£)1 

J 0 (Ro)J 



The temperature distribution in the pipe is assumed 
to have the following form: 

T(r,t) = f(r,t) + ax 



Boundary condition for the temperarure distribution is: 
T(R,t) = T w = E + ax 

f(R,t) = E 



The differential equation resulting from energy considerations 

au(r - t) 



14 1 =rf 

V b t 



The ultimate differential eouation to be solved is: 



i 

V bt 



= of 



i *f if] 

r b? ^r 2 J 



au(r , t) 



Assume the solution of the differential equation if of the 
form: 



where. 



f (r,t) = f 1 (r,t) + 
f -j (r ,t) = -lAe lnt 1 

f 2 (r,t) = -iBe int 1 



+ C 



_ £ollAll 

Jo(Htf) J 

_ lokll] 



v. r here 



f =V.-in 
1 2^ 
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The solution is: 



f ( -p + \ — «. 1 

" n|>i^2 _ n| 



_ Jn (ri )~l e int 

J 0 (R<J)J 



fo (r ,t ) = n 

n i-^ 2 - # 1 



1 - 2o(rf)l e lnt . E 
Jo(Rf)J 



The temperature distribution within the pipe is: 



T(r,t) = 



0 



i - £°i£il1 + — 
J 0 (R^)J A 



e 



1 - J o (rt) 



J 0 (Rf). 



,int 



+ E + ax , where (S = , A = - i ^ - — , © 

3* 



n ° n 

The heat transfer from the wall surface to the vapor 
in the pipe is: 



Q = - k 21Trl 



*T\ 

cTrj r = 



= R 



Therefore, the heat transfer term sought is: 



Q =-2rrklR ~ 

A 



'JJ'MI +e f JiMll e int 
. Jo(R< 5) Jo(Rf)J 



is : 



The average velocity of the vapor entering the pipe 
R 



•* / 4. \ _ Jo 2 tfru( r , t)dr 

( J " tfR 2 



v(t)=-i5 

n 



1 - 



2 Jo(rD; 

r3 j^WT 



e 



int 



Now both the velocity term and heat transfer terms 
are functions of time and can be substituted into equation 
(10) of section 11 . 24 ? 



§f = c[f# f-v(t) + Q(t) 
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Upon substitution and integration with respect to time, 
and attaching significance to the real part of the result, 
the pressure as a function of time is: 




+ ef £lllL£l 1 sin nt 

f Jo(R?) J 



The average velocity has the following form when 
returned to real notation: 



v(t) = £ 
n 



_2 Ji (RS) 

Jo im 



sin nt 



The above pressure history of the compressible volume 
and the average velocity has gouping of parameters which 
appear similar to those found in the rectangular co- 
ordinate analysis. 
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